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ABSTRACT 

We present results from the first successful millimetre (combined 33 GHz and 
35 GHz) observations of the neutron star X-ray binary Circinus X-l, using the Aus- 
tralia Telescope Compact Array. The source was clearly detected in all three observing 
epochs. We see strong evidence for a periastron flare beginning at MJD 55519.9 ± 0.04 
with estimated peak flux densities of up to 50 mJy and which proceeds to decline over 
the following four days. We directly resolve jet structures on sub-arcsecond scales. Flux 
density variability and distance from the core of nearby components suggests recent 
shock re-energisation, though we are unable to directly connect this with the observed 
flare. We suggest that, if the emission is powered by an unseen outflow, then a phase 
delay exists between flare onset and subsequent brightening of nearby components, 
with flows reaching mildly relativistic velocities. Given resolved structure positions, in 
comparison to past observations of Cir X-l, we find evidence that jet direction may 
vary with distance from the core, or the source's precession parameters have changed. 
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1 INTRODUCTION 

Circinus X-l (Cir X-l) is a confirmed (Linares et al. 2010) 
neutron star X-ray binary (NSXRB) known for its regular 
16.6 day periastron flares (radio: Whelan et al. 1977, IR: 
Glass 1978, X-ray: Tennant, Fabian & Shafer 1986), and for 
being one of the few NSXRBs to show relativistic jets, re- 
solved at a variety of scales and wavelengths (X-ray: Heinz 
et al. 2007; Soleri et al. 2009, radio: Stewart et al. 1993; 
Fender et al. 1998; Tudose et al. 2006). In the past, Cir X- 
l's flares were found to precede brightening of nearby ejecta, 
which was interpreted as re-energisation via interaction with 
unseen outflows. The time delay between radio core flar- 
ing and re-brightening of the downstream material indicated 
Lorentz factors of T > 15 (Fender et al. 2004), and while this 
has been corroborated via further advanced analysis of the 
same data-sets (Tudose et al. 2008, henceforth Tu08), sub- 
sequent observations of the source have yet to yield similar 
results (Calvelo et al. 2011 - in press, henceforth Call). 

Currently, there remains uncertainty on the orientation 
of the system and its jets. The ultra-relativistic flow veloc- 
ities would imply a jet inclination very close to the line of 
sight (6» < 5°: Fender et al. 2004). However, fitting of blue 
and red-shifted X-ray emission features (possibly originat- 
ing in the jets) have been used to calculate lower jet speeds 
(0.08c) and a jet inclination near perpendicular to the line of 



sight (Iaria et al. 2008). Furthermore, there exist multiple 
pieces of evidence for an edge-on accretion disc, including 
X-ray P-Cygni profiles (Brandt & Schulz 2000) and X-ray 
dips (Shirley, Levine, & Bradt 1999). These results need 
not be contradictory, if we assume that the jets are either 
misaligned with the orbital plane (e.g. Maccarone 2002), fol- 
low significantly non-linear flow paths (already indicated on 
larger scales in radio maps; Tudose et al. 2006), or precess. 

Tu08 reported little variation in the jet axis over a 
decade of observations (1996-2006) implying minimal pre- 
cession in the system over long periods of time. In contrast, 
more recent observations (2009/2010: Call) suggest that the 
source's 3-dimensional jet orientation may now differ from 
that observed over the 1996-2006 period (projected angle 
difference of 40 ± 5° and further from the line of sight). 
Call also found that the most significant structural vari- 
ability appeared to the north-west of the source (previously 
a region of minimal activity and speculated to be the direc- 
tion of the receding jet: Tu08), as well as structures within 
the nebula reminiscent of the filaments observed in nebula 
W50, surrounding SS 433 (Dubner et al. 1998). Call pos- 
tulate that the system's jets may have shifted recently, or 
even precess in a similar manner to those of SS 433. 

The system's quiescent and flare levels have varied sig- 
nificantly since discovery, with radio flares reaching above 
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Table 1. Cir X-l observations. The table lists the date, start Modified Julian Day (MJD) of each observation (on-source), total on- 
source time before and after flagging, MJD end time after flagging, daily measured peak flux densities (33-35 GHz) for the core and 
north-western component, and noise levels of the images from which the measurements were taken. 

Date MJD On-source Post-flag Post-flag Core Component N rms 

(UT) start time (h) on-source time (h) MJD end (mjy beam -1 ) (mjy beam -1 ) (mjy beam -1 ) 

2010 Nov 19 55519.841 8.86 6.08 55520.375 7.14 2.11 0.23 

2010 Nov 21 55521.836 8.42 3.31 55522.021 1.66 0.52 0.08 

2010 Nov 23 55523.804 8.14 5.47 55524.240 0.25 0.10 0.02 



1 Jy in the late 1970s (Haynes et al. 1978), but declining 
since 1997 to reach only 10s of mjy (Fender, Tzioumis, & 
Tudose 2005). Though inter-flare monitoring indicates Cir 
X-l remains historically 'faint', recent activity (X-ray flares 
greater than 100 mCrab: Nakajima et al. 2010, radio flares 
greater than 0.1 Jy: Calvelo et al. 2010) has renewed ob- 
servational interest in the object. Based on recent observa- 
tions with the Australia Telescope Compact Array Broad- 
band Backend (ATCA-CABB) in Call, spectral index esti- 
mates indicated Cir X-l would likely be detectable at mm 
wavelengths. NSXRBs had yet to be observed successfully 
at this frequency; thus, a detection would not only fill an un- 
explored spectral region but also allow for higher resolution 
images than had previously been achieved with ATCA. 



2 OBSERVATIONS & DATA REDUCTION 

Observations were carried out on 2010 Nov 19-20, 21-22 and 
23-24 using the ATCA-CABB in 6A configuration (mini- 
mum baseline of 337m, maximum of 5939m) at both 33 and 
35 GHz (details in Table 1). PKS B1253-055, PKS B1934- 
638, and PKS B1511-55 were used as the bandpass, flux and 
phase calibrators respectively. Between 8 and 9 hours of time 
was spent on-source each day (predicted rms noise of 26 - 
28 /xJy at both frequencies - average weather conditions). 
All data and image processing was carried out in MIRIAD 
(Sault, Teuben and Wright 1995). 

Since millimetre observations are more susceptible to 
atmospheric effects than longer wavelengths, phase stability 
had to be closely monitored via the ATCA "seeing" monitor, 
which displays rms path length fluctuations over 10 minute 
samples in microns (Middelberg et al. 2006). Phase stability 
declined as runs progressed, due to rising temperature and 
humidity. Though a switch to lower frequencies was never 
warranted (rms path length fluctuations tended to remain 
below 700 microns), once analysis began it was found that 
large segments towards the latter half of each observation 
did not have sufficient phase stability to be reliable. 

Due to phase stability decay, it was necessary to define 
the level of flagging that would improve final image fidelity. 
By dividing observations into segments, and comparing the 
resulting maps to rms path length fluctuations (VARPLT), 
we found that segments with low rms showed the expected 
point source at Cir X-l's location. Higher rms segments cor- 
responded with increasing noise towards Cir X-l's position, 
with no single distinguishable point source. In general, seg- 
ments whose majority of time was spent with rms higher 
than 400 microns failed to produce the defined structure. 
Thus for the remainder of analysis, we only use data whose 
rms values remained below 400 microns. 



3 ANALYSIS & RESULTS 

To produce maps with high signal to noise ratios, we com- 
bined 33 and 35 GHz data during the inverse Fourier trans- 
form of calibrated uv data (i.e. included both frequency's 
visibility files with multi-frequency synthesis in MIRIAD's IN- 
VERT step) . A complete spectral analysis will be covered in 
a future publication. 

3.1 Light Curves 

Declining phase stability prevented production of complete 
sets of image plane measurements for each day's full obser- 
vation. However, having produced maps for each day, it is 
apparent that Cir X-l is the only visible source, and thus 
we may estimate flux density variability from uv data am- 
plitude plots. We are, however, still unable to distinguish 
whether variation arises from the 'core' of Cir X-l (desig- 
nated as a point source situated at Cir X-l's well established 
co-ordinates: Tu08) or surrounding structure (section 3.2). 

Figure 1 shows each observation's visibility amplitudes, 
with available image plane fits over-plotted. The most no- 
ticeable feature is the multi-peaked rise in the first run's am- 
plitudes, which in comparison to the two subsequent runs, 
suggests flaring activity. A rise is also noticeable in the image 
plane fits, though it only includes a few points for each run. 
The multiple peaks may be attributed to real variations from 
Cir X-l (multi-peaked flares have been seen before: Thomas 
et al. 1978, Tu08), or resolved structure. In this case the 
radio data alone is insufficient to eliminate either possibil- 
ity. Our inability to produce high fidelity maps for the time 
periods of highest amplitude means we cannot accurately 
measure the peak flux density reached by this flare, how- 
ever, if we assume the relative difference between mean av- 
erage amplitude and core image plane values remains similar 
throughout the first observation, we may crudely estimate 
that flux densities could have reached up to 50 mjy. 

The second run's amplitudes indicate a post-flare de- 
cline, although the two image plane measurements show a 
slight rise. With this observation data set having been the 
most affected by phase disturbances, it is difficult to support 
either claim. The final observation set shows the lowest am- 
plitudes and flux densities, with the source having likely re- 
turned to quiescent levels. The flux densities measured from 
the core on each individual day are shown in Table 1 and 
also show the decline in flux across the three observations, 
supporting much of what is observed in the plots. 

3.2 Imaging 

Individual image contour plots for each day are shown in 
Figure 2. Cir X-l is resolved in all three runs, at a position 
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Figure 1. Circinus X-l ATCA-CABB rms path length fluctua- 
tions (top), ill) amplitude, and image plane fit plots (three obser- 
vations over 5 days). Measurement errors {lor) are included in the 
image plane values but are too small to appear clearly. 

that shows no shift greater than the positional errors asso- 
ciated with each set of images. Jet like structure is visible 
along a single axis (north-west to south-east) which is sim- 
ilar to that seen in previous mapping of the source at cm 
wavelengths (Tudose et al. 2006, 2008; Call). 

When interpreting the images, we must remain aware 
of the fact that we have imaged the source in both a state 
of rise and of decay, violating one of the major assumptions 
in aperture synthesis. As a result, artefacts may have been 
introduced into the image which can imitate the appearance 
of symmetrical jets. Simulations were created based on the 
visibilities from the first observation where we assumed all 
emission was from a single point source at Cir X-l's location 
with flux density S„ ~ average visibility amplitude. Images 
produced from the simulated data did indeed show signs 
of axial 'spokes' from the central source during peaks in 
amplitude; however, the intensity of these spokes beyond 
0.5 arc-seconds remained below 10 per cent of the core flux 
density (i.e. lower than the 3er noise levels measured in our 
images) and the more intense residuals left within 0.5 arc- 
seconds could be identified by their symmetrical layout and 
pairing with nearby negative regions of similar intensity. 

While we acknowledge that our first epoch image must 
be affected by the above artefacts to some degree, we re- 
main convinced that much of the resolved structure cannot 
be a result of this for a number reasons. Firstly, as per the 
simulations, the artefacts should appear symmetrical about 
the point of variation, and while there is a degree of sym- 
metry in the visible components, two distinct asymmetric 
features arise repeatedly in all images. These 'extensions' 
to the south-east and north-west of the core appear to dif- 
fer in both structure and intensity, and a separate emission 
component appears 0.75 ± 0.25 arc-seconds to the north- 
west of the core (slightly below the 'jet' axis). Secondly, the 
flagged data sets all share a similar starting hour angle, and 
though the second flagged data set is significantly shorter 
than the others, the position angle of the resulting beams 
(-30.7°, -40.9°, -41.1° for Nov 19, 21 and 23 respectively) 
remains similar across all data-sets. With this in mind, we 
would expect the axis of visible artefacts to differ between 
the first and second runs, as we are dealing with opposite 
types of light curve behaviour - rise and decay - and yet 
they remain similar. Thirdly, the final data set shows little 
or no evidence for continuing flux density decay, and as such 
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Figure 2. Radio contour maps of Circinus X-l, created using 
combined 33 and 35 GHz data. Weighting uses a robust factor of 
0.5 (an optimal compromise between natural and uniform weight- 
ing). Contours are set at -6,-4.5,-3,3,4.5,6,9,12,15,21,27,36 X rms 
noise of each epoch. Beam sizes are 0.44x0.24, 0.76x0.21, and 
0.39x0.26 arcsec for the 19th, 21st and 23rd respectively. The 
right panels have had a fitted point source subtracted at the core's 
position (the cross). The dashed circle marks component N. 



should be free of the above effects. Yet, we continue to ob- 
serve elongation of the core and the north-component in the 
final day's images. 

There is a further concern in that data whose level of 
phase errors gradually vary could produce similar artefacts 
to those from source variability. Though we have taken steps 
to eliminate afflicted data, this does not mean our images 
are not affected to some degree (see Taylor, Carilli & Per- 
ley 1999, chapter 28). The random nature of these errors 
means that over a large enough integration, artefacts should 
still appear symmetrical, making the first argument of the 
previous paragraph, coupled with multiple detections, our 
strongest test for validity; i.e. it is unlikely that persistent 
asymmetrical structure will be an artefact. 

To further examine resolved structure we performed 
point source subtraction centred at the core, with the results 
shown in Figure 2 (right hand panels). The map of the 19th 
(first epoch) shows significant residual structure, but since 
it is the day of the flare, artefacts caused by source variabil- 
ity are more likely to occur. Signs of such effects appear in 
proximity to the subtraction point as roughly symmetrical 
negative zones along an axis perpendicular to that of posi- 
tive structure. Even assuming the presence of artefacts, they 
cannot account for the brighter south-eastern core residual 
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Figure 3. Radio contour map of Circimis X-l, created using com- 
bined data from 2010 Nov 21 and 23 (33 and 35 GHz). Natural 
weighting is used for higher sensitivity at the cost of slight reduc- 
tion in beam resolution (beam size = 0.50x0.29 arcsec). Contours 
are set at -6,-4.5,-3,3,4.5,6,9,12,15,21,27,36 X rms noise of 53 fiJy. 
Layout matches that used in Figure 2. 

(approximately 50 per cent brighter; a difference of 0.8 mjy) 
nor the three other separate emission components in the 
north west quarter. Out of these remaining bright compo- 
nents, only the closest to the core (marked by the dashed 
circle in Figure 2, henceforth referred to as component N) 
remains in all subsequent images. The two structures farther 
from the core may be real and simply fade between the first 
and second observations (perhaps related to an earlier ejec- 
tion event) but without multiple detections this is difficult 
to confirm. 

The second map set also shows both south-eastern resid- 
uals and component N. We observe a streak artefact running 
the length of the image, parallel to the apparent jet axis, 
which we were unable to eliminate. This appears to be the 
source of several weak structures across the length of the 
axis, and also affects the southern extension. 

The final observation continues to show component N 
and faint extension residuals; however, the northern exten- 
sion now appears slightly offset from the 'jet' axis (further 
south than earlier examples) and is the more dominant of the 
pair of core residuals, though the difference in flux density is 
sufficiently low for noise to account for it. With variability 
likely to have reached minimum levels (as indicated by the 
image plane measurements), if any artefacts now exist in the 
map they should be a result of phase errors. This cannot be 
ruled out a cause of the low level extensions near the core. 

Figure 3 shows naturally weighted images from com- 
bined 21st and 23rd observations (post flare), resulting in a 
greater than 6<r detection of component N. Accepting com- 
ponent N as a real structure, it may be moving ejecta from 
Cir X-l or a shock caused by re-energisation of slow material 
by an unseen flow (as in Fender et al. 2004). We find that 
its flux density declines over the observations (see Table 1), 
though we cannot determine the onset of a rise as the object 
is difficult to detect on shorter data intervals. LBA observa- 
tions in our partnered submission, Miller- Jones et al. (sub- 
mitted to MNRAS, henceforth MJ2011), suggest the pres- 
ence of objects with proper motions of 35 mas day -1 . Ejecta 
moving at such speeds in our observations would result in 
detectable shifts of position between our initial and final 
images. Analysis of component N's position in our three im- 
ages shows a 0.29 ± 0.13 arc-second shift, corresponding to 
a proper motion of about 150 ± 70 mas d~ to the south be- 
tween the first and second epochs. This is the only significant 
change as all other shifts fall below respective fitting error 



limits. The shift actually brings the component closer to the 
core, perhaps implying a bent outflow or a sequence of differ- 
ent fading and brightening jet components. However, being 
a marginal difference and including the possibility that im- 
age errors have affected the component's position (i.e. more 
significant in the first day image) makes it difficult to justify 
further analysis. 

Furthermore, based on revised ephemeris estimates for 
the system (George Nicolson, private communication) our 
epoch's flare should have occurred at MJD 55519.94; about 
an hour after the rise in the amplitude plot begins. Thus, 
component N's distance from the core prohibits associa- 
tion with our flare, as it would imply a proper motion over 
2000 mas d -1 . This argument would also be true for re- 
energisation via internal shocks. Such constraints, however, 
do not apply if the event is connected to an earlier flare. 
Comparison to observations of mm flares from black hole 
(BH) XRB GRS 1915+105 (Fender & Pooley 2000) sug- 
gests that ejecta emission declines rapidly at mm wave- 
lengths (GRS 1915+105 events have duration shorter than 
1 hour), making an ejected component's possible life time of 
over 15 days questionable. Even if we assume a shallower de- 
cay, extrapolation of component N's decay rate back towards 
the earlier flare would imply minimum flux densities greater 
than 50 mjy at launch; comparable to the peak core flare lev- 
els estimated from the amplitude plots. Shocks do not suffer 
such limitations, as the energising flows may travel unde- 
tected before finally brightening on interaction with slower 
media downstream (as is observed in other BHXRBs, e.g. 
XTE J1550-564: Corbel et al. 2002 and NSXRBs, e.g. Scor- 
pius X-l: Fomalont, Geldzahler, and Bradshaw 2001). Thus, 
if we assume a single orbit's delay (16.6 days) in component 
N's variability, we estimate proper motion for a flow of ap- 
proximately 60 mas d _1 . 

One may also be tempted to carry out similar calcula- 
tions based on the slight shifts in position observed in the 
southern and northern extensions, but the possibility that 
they are affected by any of the artefacts we have discussed 
would severely limit our confidence in such an exercise. 



4 DISCUSSION 

Having established that image artefacts predominantly af- 
fect the structures nearest to the core, our best jet axis 
estimate uses component N's position, giving an angle of 
121 ± 5° E of N. Though this estimate falls within the 
range quoted in Tu08 (129 ± 13° E of N) and is only 
marginally steeper than the estimate from milli-arcsecond 
structure in MJ2011's LBA observations (113° E of N), it is 
approximately 50° shallower than the axis estimated from 
5.5 GHz residual emission near the core in Call. Further- 
more, whereas MJ2011 observe symmetrical jets, our maps 
indicate a scenario much like that observed in Call and 
Tu08 where asymmetry of near-core residuals indicates a 
brighter southern jet. Like Call, we also observe distinct 
activity from a region that has previously been associated 
with the receding jet (Tu08) in the form of component N. 

The fact that the various structural orientations appear 
on different scales across observations taken within months 
of each other supports a scenario where the jet angle varies 
as we move further from the system's core. Starting off at 
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relatively high angles to the line of sight and near E-W as 
indicated in MJ2011, but steadily turning N-S and closer 
to the line of sight as we move down the flow to arcsec- 
ond scales as in Call. The origin of such deviation may be 
precession, or quasi-static jet kinks. The unexpected appear- 
ance of structures to the north west (component N in this 
case) may be a result of asymmetry in the density of inter- 
acting media, or deviations in a jet close to the line of sight 
(as suggested in Call). 

The jet velocities reported in MJ2011 are of similar or- 
der to our estimates for energising flows in section 3 (when 
we assume a full orbit's phase delay), both of which fall 
far below the ultra-relativistic speeds calculated by Fender 
et al. (2004). However, we too cannot completely rule out 
the existence of such ultra-relativistic flows. There remains 
the possibility that we did not observe the full flare event 
associated with our specific orbit (ephemeris estimates are 
not perfect, nor every flare event the same), altering the ref- 
erence time on which our component N calculation is based. 

Reconciling these combined results with historic ob- 
servations such as the ultra-relativistic velocities of Fender 
et al. (2004) and lack of precession observed in Tu08 is diffi- 
cult. Lacking details on sub-arcsecond structure at the time 
of the Tu08 observations, we cannot establish whether Cir 
X-l's jets have always exhibited variable inclination with 
distance from the core (though jet curvature in large scale 
images from Tudose et al. 2006 indicate some level of this) . 
By grouping radio observations made after 2008 and com- 
paring them to older results, we begin to observe an almost 
bimodal divide in Cir X-l's behaviour, with lower jet ve- 
locities, precession (indicated by the difference in jet angle 
observed at 5.5 GHz: Call) and multiple incidences of activ- 
ity to the north-west of the core rather than the south-east 
being recent changes. This may indicate that, along with 
the decline in intensity over the last decade, the system's jet 
structure can also vary significantly with time. 

Such a scenario draws further similarities between 
Cir X-l and SS433/W50. Hydrodynamical simulations by 
Goodall et al. (2011) show that the structures visible in neb- 
ula W50 cannot be reproduced by extrapolation of SS433's 
current jet activity back through time. Instead they show 
that at least three distinct epochs of jet activity, each with 
a different level of precession, were required to produce the 
nebula's layout. If the filaments visible in Call's Cir X- 
1 nebula maps are indeed caused by jet-ISM interactions, 
then the precession angle of the jets would have had to be 
significantly larger than that indicated by pre-2008 radio ob- 
servations. Thus, variability of Cir X-l's jet characteristics 
should not be entirely unexpected. 



5 SUMMARY 

We have detected Circinus X-l and sub-arcsecond structure 
at millimetre wavelengths for the first time. Light curves in- 
dicate the rise phase of a periastron flare event occurred on 
Nov 19, decaying over the following four days. Radio maps 
show resolved structure near the core and a 6<r level com- 
ponent to the north-west whose position angle, though in 
agreement with the jet axes of Tu08 and MJ2011, indicates 
activity in a direction previously associated with the reced- 
ing jet. Variability of the north-western component, likely 



the result of re-energisation by unseen outflows, cannot be 
physically reconciled with the observed flare but if caused by 
an immediately previous periastron event implies a proper 
motion of 60 mas day -1 . Given the observed differences in 
jet behaviour compared to results prior to 2008, we believe 
that this work, in conjunction with results from Call and 
our partnered submission, MJ2011, suggest Cir X-l's jet be- 
haviour; including flow velocity and precession parameters, 
can vary on secular timescales. 
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